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FIRST ORDER SYSTEM 
(floating point) 



PURPOSE: 



To solve the equation, 

c - r = T (dr/dt) 
for r, where c is given as a function of time. 
METHOD ; 

The following recursion relation is used: 

r = A(c -c ,) + B(r , -c ,) + c 
n n n-1 n-1 n-l n 

where: A = -~ (1- c ' ), B = e 

z^t ' 

giving the value of r at the n'th instant in terms of c at the n'th instant and 
the previous r and c at the (n-l)'th instant. Since all storages are upgraded 
internally, it is only necessary to enter the routine with successive c's, 
exiting with the successive r's, if the interval between n and n-1 is unchanged. 
The routine assumes that c varies linearly with time between points, so 
that it is only necessary to choose the points on c so that c is sufficiently 
well represented as a series of straight lines. 

RESTRICTIONS: 



CALLING SEQUENCE: 

TRA (location of this routine); exit is to next half-word. On entry, 

A and R should contain c ; on exit, A and R will contain r . 

n n 



4 A " f w ? f 
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PARAMETER STORAGE ; 

The following locations must be filled as indicated with the correct 
initial conditions before this subroutine is used. The r and c values are 
upgraded automatically; time is assumed to increase by an amount z^t 
each time through the routine. 

Lo 4 30 r Lo + 34 At 

Lo + 32 c Lo + 36 T 

In addition, the following are calculated and made available in the locations 
shown, but need not be filled with initial conditions. 

Lo + 42 c-r Lo +44 dr/dt 

New r or c . values can be stored in the indicated locations at any time, 

n-X n-i 

resulting in a discontinuity; thus discontinuities in either r or c can be 
solved for accurately. 

Z^t may be changed as desired to provide larger or smaller steps in 
time, to show more or less detail, or to suit the nature of the input function 
c. T may also be changed in a step-wise fashion, thus yielding an approxi- 
miate solution to many non -linear equations. 

ROUTINE REQUIRED : 

Exponential (AN -044) in 1050; called for by TRA instruction in (Location 
of this routine + 46) 

STORAGE REQUIRED; 

100 _. sectors, plus use of both high speed loops, 
(o) 

REMARKS ; 

The solution is "exact" in that the recursion relationship will give the 
exact solution (except for round-off errors) if c varies linearly between 
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points and tIq a constant. Non-linear and tixne-varying solutions can be 
approximated by choosing At small enough so that ^ is essentially constant 

over each computation Interval. The constants A and B in the recursion 

-10 
relationships are computed with an accuracy within 0. 5 x 10 



TIME: 



0. 8 seconds to solve for r and evaluate A and B. 

n 

0.2 seconds to solve for r (no change in £^t on ). 

n 



CHECKOUT: 



This routine generated correct solutions to: 

1 -tA 

(1) simple step response: r = -7 — -rr— = 1 - e 

*^ ^ *^ 3(s + l) 

(2) non-linear system step response: r = 



s(l+8+t8) 1+t 

FI LE UNITERMS ; 

First Order System, Differential Equation, Time Constant 
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